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A  new  fluid  theory  in  the  guiding-center  and  gyrotropic  approximation  derivable  from  the 
ensemble-averaged  Vlasov-Maxwell  equations  that  included  the  effect  of  wave-particle  interactions 
for  weakly  turbulent,  weakly  inhomogeneous,  nonuniformly  magnetized  plasma  was  recently  given 
by  Jasperse,  Basu,  Lund,  and  Bouhram  [Phys.  Plasmas  13,  072903  (2006)].  In  that  theory,  the 
particles  are  transported  in  one  spatial  dimension  (the  distance  s  along  the  magnetic  field)  but  the 
turbulence  is  two-dimensional.  In  this  paper,  which  is  intended  as  a  sequel,  the  above  theory  is  used 
for  quasisteady  conditions  to  find:  (1)  a  new  formula  for  the  perpendicular  ion-heating  rate  per  unit 
volume  VV/±0s;),  where  Wi±(s)  decreases  for  large  5  by  what  we  call  the  “finite  ion  gyroradius 
effect”;  and  (2)  a  new  formula  for  the  perpendicular  ion  temperature  at  low  altitudes,  Ti±(s). 
Parametrized  calculations  for  Ti±(s)  are  also  given.  ©  2006  American  Institute  of  Physics. 

[DOI:  10.1063/1.2364475] 


I.  INTRODUCTION 

In  Ref.  1,  hereafter  referred  to  as  I,  a  new  fluid  theory 
was  presented  in  the  guiding-center  and  gyrotropic  approxi¬ 
mation  that  included  the  effect  of  wave-particle  interactions 
for  weakly  turbulent,  weakly  inhomogeneous,  nonuniformly 
magnetized  plasma.  It  is  the  object  of  the  present  paper  to 
use  the  new  theory  for  quasisteady  conditions  to  derive  a 
formula  for  the  perpendicular  ion-heating  rate  per  unit  vol¬ 
ume  WiL(s),  and  the  perpendicular  ion-temperature  profile  at 
low  altitudes  Ti±(s),  for  a  downward  Birkeland  current  re¬ 
gion.  The  present  paper  is  intended  as  a  sequel  to  I. 

The  Birkeland  currents  are  a  system  of  upward  and 
downward  magnetic  field-aligned  electrical  currents  in  the 
auroral  zone  that  flow  between  the  magnetosphere  and  the 
ionosphere  where  particle-number  flux,  momentum  flux,  and 
energy  flux  are  exchanged.  For  an  early  study  of  the  statisti¬ 
cal  properties  of  the  Birkeland  currents  using  TRIAD  satel¬ 
lite  data  (see  Ref.  2). 

Over  the  years,  much  progress  has  been  made  in  iono¬ 
spheric  plasma  physics  (see,  for  example,  the  recent  text  by 
Schunk  and  Nagy  '  and  the  references  therein).  The  plasma  in 
most  of  the  Earth’s  ionosphere  is  collisional  and,  as  a  result, 
sometimes  not  too  far  from  equilibrium.  By  the  same  token, 
the  use  of  magnetohydrodynamic  (MHD)  theory  has  greatly 
increased  our  understanding  of  solar  wind  and  magneto- 
spheric  plasmas,  which  are  essentially  collisionless.  (See  the 
text  by  Parker4  and  the  references  therein  for  a  discussion  of 
the  progress  that  has  been  made  in  this  area.)  However,  it  is 
well  known  that  when  an  E  field  parallel  to  B  develops  in  the 
Earth’s  magnetosphere,  MHD  theory  is  no  longer 
applicable  4  Furthermore,  magnetospheric  plasma  is  far  from 
equilibrium  and  is  almost  always  highly  turbulent.  Since 
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MHD  theory  cannot  be  used  to  treat  the  field-aligned 
Birkeland  current  system,  a  more  general  approach  is 
needed. 

In  I,  two  traditional  approaches  to  plasma  turbulence 
were  briefly  discussed.  For  a  review  of  and  references  to  the 
theoretical  literature  in  this  area,  which  will  not  be  repeated 
here  (see  the  Introduction  in  Ref.  1).  There  is  a  vast  literature 
on  upward  auroral-current  regions,  which  also  will  not  be 
discussed  here  (see  Ref.  5  for  a  review  of  that  literature). 

The  properties  of  downward  auroral -current  regions  are 
less  explored.  They  have  been  observed  by  the  ISIS-2,6 
DE-l,7  S3-3,8  FREJA,9  FAST,10"13  CLUSTER,14  and  other 
satellites.  From  these  studies,  we  see  that  downward  auroral- 
current  regions  are  often  characterized  by 

(1)  diverging  electrostatic  shocks,  downward  pointing 

(2)  upflowing  low-energy  (<few  keV)  field-aligned  elec¬ 
trons; 

(3)  intense  transverse  ion  heating,  energetic  ion  conics; 

(4)  small-scale  density  cavities; 

(5)  intense  extremely  low  frequency  and  very  low  frequency 
electric-field  turbulence  and  fast  bipolar  solitary  struc¬ 
tures; 

(6)  small  fluxes  of  energetic  (>few  keV)  magnetospheric 
particles. 

Some  theoretical  work  has  been  carried  out  using  the  fluid 
equations  for  the  current-driven  ion-cyclotron  instability, 1  ^ 
particle-in-cell  simulations  in  the  electrostatic 
approximation,16-20  and  velocity-space  diffusion  models  for 
the  problem.21-24 

In  Refs.  21-24,  for  lack  of  detailed  satellite  measure¬ 
ments,  it  was  assumed  that  magnetospheric  ions  were  heated 
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in  downward  current  sheets  by  cyclotron  resonance  with 
electromagnetic  turbulence  at  the  ion  cyclotron  frequency 
ft,,  Now  that  we  have  detailed  satellite  measurements,  we 
are  able  to  determine  that,  at  least  for  the  three  cases  of 
FAST  satellite  data  studied  in  this  paper,  the  turbulence  near 
ft,  and  its  harmonics  is  electrostatic.  Since  the  parallel  wave- 
number  for  electrostatic  turbulence  is  much  larger  than  that 
for  electromagnetic  turbulence,  an  integration  over  the  paral¬ 
lel  ion  velocity  must  be  performed  in  order  to  calculate  the 
turbulent,  perpendicular  ion-heating  rate  per  unit  volume, 
Wi±.  This  complicates  the  formula  for  WiL  (see  Sec.VI). 

The  question  arises:  Is  the  one-dimensional  fluid  theory 
reviewed  in  Sec.  II  of  this  paper  and  derived  in  I  applicable 
to  downward  Birkeland  current  regions?  The  answer  is  a 
qualified  “yes.”  This  issue  is  discussed  in  detail  in  Sec.  Ill, 
but  it  suffices  to  say  here  that,  as  long  as  we  apply  the  one¬ 
dimensional  theory  to  downward  current  sheets  away  from 
the  edges  of  the  small-scale,  “U-shaped,”  quasistatic  poten¬ 
tial  structures  where  the  EiXB  drift,  udJ  is  small  compared 
to  the  ion  drift  along  B,  uh  the  one-dimensional  theory  ap¬ 
plies.  Near  the  quasistatic  potential  structures,  m,  and  ud  are 
comparable,  and  the  one-dimensional  theory  does  not  apply. 
Here,  a  two-dimensional  theory  must  be  used. 

As  discussed  in  I,  if  we  could  find  a  renormalized  solu¬ 
tion  of  the  wave-kinetic  and  plasma-kinetic  equations,  the 
problem  would  be  solved.  The  renormalized  propagator 
could  then  be  used  to  find  the  renormalized  distribution  func¬ 
tions,  dielectric  screening  function,  and  the  spectral  density 
for  the  longitudinal  electric-field  fluctuations  for  the  turbu¬ 
lent  plasma.  However,  the  development  of  a  renormalized 
solution  for  turbulent,  inhomogeneous,  nonuniformly  magne¬ 
tized  plasma  is  indeed  a  formidable  problem.  The  idea  that 
we  present  in  this  paper  is  to  bypass  this  difficult  problem  by 
solving  the  fluid  equations,  where  satellite  measurements  can 
be  used  to  specify  the  turbulent,  electric-field  fluctuations 
and  where  satellite  measurements  of  the  low-order  velocity 
moments  of  the  particle  distributions  can  also  be  used  to 
specify  the  initial  and  boundary  conditions. 


II.  GYROTROPIC  GUIDING-CENTER  FLUID  THEORY 
FOR  THE  BIRKELAND  CURRENT  SYSTEM 

In  I,  we  gave  a  new  fluid  theory  in  the  guiding-center 
and  gyrotropic  approximation  derivable  from  the  ensemble- 
averaged  Vlasov-Maxwell  equations  that  included  the  effect 
of  wave-particle  interactions  for  weakly  turbulent,  weakly 
inhomogeneous,  nonuniformly  magnetized  plasma.  It  was 
assumed  that  the  turbulence  is  random  and  electrostatic,  and 
that  the  Fokker-Planck  operator  could  be  used  to  calculate 
the  correlation  functions  that  describe  the  turbulent,  wave- 
particle  interactions.  At  present,  the  theory  has  been  worked 
out  in  one  spatial  dimension,  i.e.,  the  distance  s  along  the 
geomagnetic  field,  but  the  turbulence  is  two-dimensional 
since  electric-field  fluctuations  that  are  both  perpendicular 
and  parallel  to  the  geomagnetic  field  are  considered.  In  this 
section,  we  review  the  theory  for  the  reader’s  convenience. 


From  I,  the  fluid  equations  in  the  guiding-center  and 
gyrotropic  approximation  for  continuity,  parallel  momentum 
balance,  parallel  energy  balance,  and  perpendicular  energy 
balance,  together  with  Poisson’s  equation,  are 


d  d 

—na+B—(nauJB)  =  0, 
dt  ds 


(i) 


d  d 

—mjiaua  +  —na(TM  +  mau2a) 
dt  ds 


1  dB  2  .  d(f> 

~~  D  »  all  +  ™ oMa  ~  ^a_L.)  Q aflat  ~  —  ^ a\\>  (^) 

B  as  ds 


1  d  .  ,  d  1  dB 

a II  +  mcXa)  +  ~  ”a4ail  "  "aWall  "  <7a±) 

2  dt  ds  B  ds 


d(f> 

+  q</icMa—  =  W^n, 
ds 


(3) 


d 

—nj, 
dt 


<2* 
ds ' 


a±  +  B-—(naqal/B-)  =  Wa±, 


(4) 


—  =  -4t t2j 


ds 


(5) 


where  s  is  the  distance  along  the  geomagnetic  field,  denoted 
by  B,  and  t  is  the  time.  Gaussian  units  are  used.  For  particles 
of  type  a,  na  is  the  number  density,  ua  is  the  parallel  drift 
velocity,  Tall  is  the  parallel  temperature,  Tal  is  the  perpen¬ 
dicular  temperature,  is  the  total  parallel  energy  flux  per 

particle,  and  qal  is  the  total  perpendicular  energy  flux  per 
particle.  By  total,  we  mean  the  sum  of  the  drift  and  random 
parts.  We  note  that  naq^(nc/iaL)  is  the  total  parallel  (perpen¬ 
dicular)  energy  flux  for  particles  of  type  a.  The  definitions 
for  q^ |  and  qal  should  not  be  confused  with  the  heat  fluxes 
per  particle,  also  denoted  by  q^  and  q±  in  Braginskii,  ^  but 
defined  differently  here.  The  electrostatic  potential,  denoted 
by  </>,  and  the  low-order  velocity  moments  of  the  one-particle 
distribution  functions  are  functions  of  5  and  r.  We  have  also 
introduced  the  quantities  and  Wal ,  where  is 

the  rate  of  transfer  of  momentum  per  unit  volume  for  par¬ 
ticles  of  type  a  due  to  wave-particle  interactions,  and 
W^iWai)  is  the  rate  of  transfer  of  parallel  (perpendicular) 
energy  per  unit  volume  for  particles  of  type  a  due  to  wave- 
particle  interactions.  The  transfer  rates  are  also  functions  of  s 
and  t.  The  are  related  to  the  anomalous  (turbulent)  re¬ 
sistivity  for  the  problem,  and  and  Wal  are  the  anoma¬ 
lous  (turbulent)  parallel  and  perpendicular  heating  or  cooling 
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rates  per  unit  volume  for  particles  of  type  a.  The  other  quan¬ 
tities  in  ( 1 )— (5)  have  their  usual  meaning.  The  wave-particle 
transfer  rates  per  unit  volume  were  defined  in  I  as 


0 


mcP  II 

mav]l2 

mav\l2 


Cn. 


(6) 


Here,  Ca  is  the  gyrotropic  average  of  Ca,  which  denotes  the 
wave-particle  correlation  function  for  particles  of  type  a. 

In  the  guiding-center  and  gyrotropic  approximation,  we 
considered  a  Cartesian  velocity-space  coordinate  system  that 
slowly  changes  its  orientation  as  s  varies  along  B  so  that  the 
vz  axis  is  always  parallel  or  anti-parallel  to  B.  We  then  as¬ 
sumed  that  the  fluctuating  electric  field  is  random  (Markov¬ 
ian),  that  the  length  and  time  scales  for  the  one-particle  dis¬ 
tribution  functions  and  the  two-particle  correlation  functions 
separate,  and  that  Ca  is  given  by  a  velocity-space  Fokker- 
Planck  operator. 26,27 


<9  ,  -  1  d  d  I 

■  -  •  (F{  +  F ?a)fa  +  . 

d\  2  dX  dx 


(7) 


where  is  the  gyrophase  angle  in  the  guiding-center  coordi¬ 
nate  system.  Here,  F{,  F£,  and  are  functionals  of  e  and 
(\SEI 2\),  where  e  and  (|<$E2|)  are  the  dielectric  screening 
function  and  spectral  density  of  the  longitudinal  electric-field 
fluctuations  for  the  turbulent  plasma,  respectively.  In  (7),  fa 
is  the  gyrophase  average  of  the  one-particle  distribution 
function  for  particles  of  type  a,  Ffa  denotes  the  friction  due  to 
the  correlated  effect  of  the  electric-field  fluctuations  on  the 
particle  orbit,  Fpa  denotes  the  friction  due  to  the  polarization 
of  the  turbulent  plasma  by  a  test  particle,  and  denotes  the 
diffusion  tensor  due  to  the  correlated  effect  of  the  electric- 
field  fluctuations  on  the  particle  orbit.  In  I,  we  also  gave 
explicit  expressions  for  Ma)|,Wa)i>  and  Wal  for  weak,  ran¬ 
dom,  electrostatic  turbulence  in  the  guiding-center  and  gyro¬ 
tropic  approximation.  They  are: 


£||  \.  -T.  -  V  d  „  .  w  /  &||  1  2/ 


+  (  ]jrn(^)4'WoIm[fi(5,r,k,(i>)  ]  f  +  A:||L?||  —  0>) , 


(8) 
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Wa±(s,t)  =  (^J  f  d3vfa(s,t,u±,v h)  J  J  ^w|(^)n[72.1(4)-y2+,(^)]<|^E2|(i,r,k,(o» 
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(10) 


where  €a=k1v±/Fla  and  VL^qJUm.f  to  include  the  sign 
of  qa.  In  addition,  Jn(g)  is  the  usual  Bessel  function  of  order 
h,  k±(k\)  is  the  perpendicular  (parallel)  part  of  k,  and 
k2-k\  +k2.  Here  k  is  a  vector  in  the  guiding-center  coordi¬ 
nate  system.  The  other  symbols  have  their  usual  meaning. 
We  note  here  that  the  following  symmetries  hold: 


I 

(|^£2|(5,r;k,cu))=(|<5£:2|(^,r;-k,-cu));  and  e(s>t\  k,cu)* 
=  e(s,f;-k,-a>).  We  also  note  that  since  B  is  a  function  of  s , 
so  are  Cta. 

Equations  (l)-(5)  are  a  set  of  4/V+  1  equations  for  the 
6N+  1  unknowns:  na ,  ua ,  Tal>  qal],  qaL,  and  </>,  where  N 

is  the  number  of  plasma  constituents.  For  the  set,  2N  closure 
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conditions  are  needed.  When  the  model  for  B  is  given,  and 
when  fa>(\SE2\)  and  e  are  specified,  then  the  set  of  fluid 
equations,  together  with  (8)— ( 1 0)  and  the  closure  conditions, 
may  be  solved  subject  to  appropriate  initial  and  boundary 
conditions  to  give  the  spatial  and  temporal  evolution  of  the 
turbulent  plasma.  Examples  of  how  this  is  done  for  quiescent 
plasma  are  given  in  Secs.  VII  and  VIII  of  I. 


III.  CONDITIONS  FOR  THE  APPLICABILITY 
OF  THE  ONE-DIMENSIONAL  FLUID  THEORY 
TO  THE  BIRKELAND  CURRENT  REGIONS 

The  question  arises:  Is  the  one-dimensional  fluid  theory 
reviewed  in  Sec.  II  and  derived  in  I  applicable  to  downward 
Birkeland  current  regions?  The  answer  is  a  qualified  “yes.” 
From  I,  we  may  identify  seven  conditions  that  must  be  sat¬ 
isfied  in  order  to  apply  the  one-dimensional  theory  to  down¬ 
ward  current  sheets.  They  are:  (1)  the  length  scales  for  the 
variation  of  B  and  E  must  be  large;  (2)  the  time  scale  for  the 
variation  of  E  must  also  be  large;  (3)  the  gyrophase  average 
approximation  must  be  valid;  (4)  the  E1XB  drifts  must  be 
small;  (5)  the  turbulence  must  be  electrostatic;  (6)  the  turbu¬ 
lence  must  also  be  weak;  and  (7)  the  Fokker-Planck  separa¬ 
tion  of  length  and  time  scales  must  be  valid.  In  this  section, 
we  discuss  each  of  these  items  in  consecutive  order. 

From  I,  we  see  that  the  use  of  guiding-center  coordinates 
implies  the  following  ordering  of  length  scales: 

hw  ^  ai\\i  Ib±  ^  ^£11  ^  ai\\\  l El  ^  ail  >  (11) 

where  the  lengths  and  ai±  are  defined  as  and 

atl  =  |ui±/n,|,  and  /  denotes  the  length  scale  for  the  B  and  E 
variations  given  in  the  parallel  and  perpendicular  directions 
by  lBx=\(B~ldB/dx)~l\  and  lEx=\(E~ldE/dx)~l\,  where  jc  is 
the  distance  either  parallel  or  perpendicular  to  the  field.  In 
the  three  FAST  satellite  data  sets  that  we  have  examined  in 
this  paper,  we  found  that  protons  dominate  the  heavier  ion 
constituents.  At  FAST  altitudes  for  the  three  data  sets,  we 
find  that  30-70  km/s,  Ti±  ~  100-200  eV,  and  B 
—  0.1  G  to  give  the  following  average  values:  a^- 40  m  and 
ai±=  1.2 X  102  m.  In  addition,  /?($)  —  B0(a/s)3,  <3  —  6.4 
X  103  km,  s2~  10.4X  103  km,  where  s2  is  the  satellite  alti¬ 
tude,  so  we  estimate  that  Igw^iw  and  lBJalL  are  both  ^104. 
Now  consider  the  E-field  length  scales.  In  general,  there  are 
two  cases:  the  large-scale  E±  that  leads  to  convection  E± 
XB  drifts  as  in  the  Polar  Cap  and  dayside  Cusp-Cleft  re¬ 
gions,  and  the  small-scale  E±  near  the  edges  and  sometimes 
in  the  middle  of  the  U-shaped,  quasistatic  potential  structures 
that  are  known  to  occur  in  both  upward  and  downward  cur- 
rent  sheets.  ’  ’  Consider  the  small-scale  E±.  There  are  two 
situations:  near  the  shock  edges  and  away  from  the  shock 
edges.  From  FAST  data  for  the  downward  current  sheet  be¬ 
ginning  at  UT  1997-01-23/14:37:40  and  extending  for  about 
10  s,  we  find  that  E^^^O.S  V/m  in  >2.5  km,  so  that 

4hrk/a,i  s  42.  (12) 

From  an  examination  of  the  time  series  of  the  fluctuating 
electric  field  along  the  satellite  track  in  a  downward  current 
sheet  for  the  time  period  from  :07:30  to  :07:40  of  orbit  1626 


shown  in  Fig.  1  of  Ref.  10,  we  see  that  for  the  no-shock 
region  £'"oshock <  1 1  mV/m  in  >50  km,  so  that 

lf[hock/ai±  >  8.3  X  102.  (13) 

Note  that  the  waveform  labeled  as  the  dc  field  in  Ref.  10, 
Fig.  1,  is  actually  sensitive  to  frequencies  up  to  1  kHz/0 
Consider  lm.  From  Ref.  24,  Fig.  6  in  the  auroral  zone, 

^noshock  js  _o  g  mV/m  jn  9X  ]02  km  SQ 

/™shock/ai||  ~  2.3  X  104.  (14) 

The  conclusion  here  is  that  the  guiding-center  coordinate 
system  may  be  used  even  near  a  small-scale,  electrostatic 
potential  structure  providing  that  the  length  scale  for  E±  is 
sufficiently  large  compared  to  <3,^. 

In  the  one-dimensional  fluid  theory  given  in  I,  we  have 
assumed  that  the  ensemble-averaged  E  (the  background  E± 
and  E,()  varies  on  a  time  scale  that  is  sufficiently  large.  From 
I,  we  see  that 

r»|2w/a,|~5X  10-3s,  (15) 

where  r  is  the  time  scale  for  the  variation  of  EL  and  E;. 
Thus,  the  ensemble-averaged  E  field  must  vary  on  a  time 
scale  that  is  large  compared  to  5  X  10"3  s. 

In  deriving  the  fluid  equations,  we  have  used  the  gy¬ 
rophase  average  of  the  ion  and  electron,  one-particle  distri¬ 
bution  functions.  In  Ref.  23,  Fig.  2,  a  gyrophase-averaged 
kinetic  calculation  for  the  ion  conic  is  compared  to  the  ex¬ 
perimental  observation  in  the  dayside  Cusp-Cleft  region  near 
6  X  103  km  from  the  Earth’s  surface  for  a  downward  current 
sheet.  We  see  that  the  gyrophase-averaged  contours  at  108 
and  107  s3/km6  are  within  about  10%-20%  of  the  experi¬ 
mental  contours. 

In  our  derivation  of  the  fluid  equations  in  one  spatial 
dimension,  we  have  neglected  the  ExXB  drift  compared  to 
the  ion  drift  along  B.  We  note  here  that  gradient  B  and  cur¬ 
vature  B  drifts  are  much  smaller  than  ExXB  drifts  for  Bir¬ 
keland  current  sheets.  First,  we  consider  FAST  altitudes. 
Again,  there  are  two  situations:  away  from  a  small-scale 
electrostatic  shock  and  near  a  small-scale  electrostatic  shock. 
From  an  examination  of  the  time  series  for  the  fluctuating 
electric  field  along  the  satellite  track  beginning  at  UT  1997- 
01-23/14:37:40  and  extending  for  about  10  s,  we  find  that 
E noshock <  j  j  my/m  jn  >50  km  and  w,^50  km/s,  so  that 

|w;/w/oshock>45,  (16) 

where  ud  is  the  E±  X  B  drift  velocity.  We  see  that  the  neglect 
of  ud  is  a  reasonable  approximation  for  the  above  data  set 
and  the  one-dimensional  theory  applies  away  from  the  small- 
scale,  quasistatic  potential  structure.  From  Ref.  10,  Fig.  2, 
orbit  1797,  we  see  that  Ef0** <  1  V/m  in  >15  km,  so 

\ui/ud\shock-l/2.  (17) 

Thus,  near  a  small-scale,  quasistatic  shock  in  the  long-range 
potential  region  of  a  downward  current  sheet  as  shown  in 
Ref.  10,  Fig.  2,  the  neglect  of  ud  compared  to  m,  breaks 
down,  and  a  two-dimensional  fluid  theory  must  be  used.  At 
CLUSTER  altitudes  for  a  downward  current  sheet,  we  see 
from  Ref.  14,  Fig.  3,  that 
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|M,/M/oshock  >  150,  (18) 

and  that 

|«,/«/hock~2.6.  (19) 


The  neglect  of  ud  for  a  no-shock  situation  appears  valid  for 
altitudes  up  to  4.3 RE  (geocentric),  but  near  the  shock,  the 
one-dimensional  theory  still  cannot  be  used  at  CLUSTER 
altitudes  for  downward  current  sheets. 

We  have  also  assumed  that  the  turbulence  is  electro¬ 
static.  In  Secs.  IV  and  VI  of  this  paper,  we  will  subtract  the 
Doppler-shifted,  quasistatic  turbulence  from  the  real,  time- 
dependent  turbulence  near  the  proton  gyrofrequency  and  its 
harmonics  and  consider  the  latter  signal  to  be  the  one  that 
heats  the  ions  by  cyclotron  resonance.  From  the  wave  data 
for  <5E  and  SB  near  H,,  we  find  that  the  Poynting  flux  enter¬ 
ing  the  flux  tube  at  s2  is  ^0.01  /zW/m2.  From  the  integrated 
total  energy  density  conservation  equation  for  steady-state 
conditions  given  by  Eq.  (41)  of  I  for  electrostatic  turbulence, 
we  see  that 


2  na{qa I,  +  qal) 

e,i 


P2  1  dB^ 


(20) 


On  the  left-hand  side  of  (20),  we  can  show  that  neqe\\  domi¬ 
nates  the  ion  total  energy  flux.  Using  the  satellite  particle 
data  at  s2  that  we  will  analyze  in  Secs.  IV  and  VI,  and  the 
approximate  theory  given  in  Ref.  24  to  estimate  the  second 
term  in  (20),  we  may  show  that  the  two  terms  on  the 
left-hand  side  of  (20)  add  to  give  —14  ytxW/m2.  Since  the 
Poynting  flux  entering  the  flux  tube  at  s2  is  very  small  com¬ 
pared  to  the  sum  of  the  two  terms  on  the  left-hand  side  of 
(20),  we  conclude  that  the  energy  flux  to  drive  the  two  terms 
on  the  left-hand  side  of  (20)  is  not  supplied  by  the  Poynting 
flux  near  d,.  As  a  check  on  this  result,  we  may  also  use  the 
satellite  data  for  j\\(s)  and  the  approximate  theory  given  in 
Ref.  24  to  estimate  El{(s)  on  the  right-hand  side  of  (20).  We 
obtain,  to  within  experimental  error, 


ds  _/||£n  =  14  ynW/rrr. 


(21) 


Thus,  we  conclude  that  the  turbulence  at  frequencies  near  d, 
and  its  harmonics  is  essentially  electrostatic,  and  all  of  the 
particle  energy  flux  to  within  experimental  error  is  supplied 
by  the  integral  of  on  the  flux  tube  and  not  by  electro¬ 
magnetic  waves  near  d,  propagating  in  from  outside  the  flux 
tube. 

In  deriving  expressions  for  the  turbulent  wave-particle 
transfer  rates  in  I,  we  have  made  the  weak  turbulence  ap¬ 
proximation.  This  implies  that  the  fluctuating  quantities  are 
ordered  by 

WJfa\  ~ \SnJna\  ~  \qaS<f>/Ta\  ~  (2m/k\D) 

X(8E2/8miaTa)ul  <g.  1,  (22) 


where  \D  is  the  Debye  length,  and  k  is  an  average  or  “char¬ 
acteristic”  wavenumber  for  the  turbulence.  Using  the  nota¬ 


tion  W=SE2/8imT,  we  see  that  for  weak  turbulence, 
W]/2<^  1.  For  the  FAST  data  set  beginning  at  UT  1997-01- 
23/14:37:40,  we  see  that  /i,  — 1  —  3  cm-3,  Ti±  —  100-200  eV, 
<5£2  — 2.1  X5.1  X  10-5V2/m2.  We  find  that  W-9.9X10-6, 
and  therefore 

W1/2  —  3.1  X  10"3<  1.  (23) 

In  (22),  k\D  is  an  unknown  parameter  that  depends  on  the 
saturated  state  of  the  turbulent,  inhomogeneous,  nonuni- 
formly  magnetized  plasma.  If  we  accept  results  from  a  par¬ 
ticle  simulation  for  simplified  boundary  conditions  in  a  ho¬ 
mogeneous  medium,18  we  can  then  estimate  that  —  k±  and 
k±\D~~0. 1,  so  that  2 1/2/&\D—  14.  As  long  as  2 U2/k\D  is  not 
too  much  larger  than  this  for  realistic  conditions  on  down¬ 
ward  current  sheets,  the  above  data  set  satisfies  the  ordering 
given  by  (22). 

The  use  of  the  Fokker-Planck  method  to  evaluate  the 
interaction  between  the  particles  and  the  turbulence  also  as¬ 
sumes  that  the  length  and  time  scales  between  the  one- 
particle  distribution  functions  (and  their  moments)  and  the 
correlation  functions  for  the  electric-field  fluctuations  sepa¬ 
rate.  In  Sec.  VI  of  I,  we  assumed  that  and 

r1>r2  — rac.  Here,  rac  is  the  autocorrelation  time  for  the 
random  electric-field  fluctuations  in  the  magnetized  plasma. 
At  FAST  altitudes,  \D<>  —  20-50  m  and  T*ac 2 tt/ O,  —  5 
X  10-3  s,  so  that  li^> 20-50  m  and  r^SXlO^s.  This 
means  that  the  low-order  velocity  moments  of  fa  may  vary 
on  a  time  scale  on  the  order  of  5  X  10-3  s  or  slower.  These 
assumptions  imply  that  if  large-amplitude  coherent  structures 
(waves  or  solitons)  at  frequencies  near  D,  and  its  first  few 
harmonics  are  present  and  r2  becomes  large  compared  to  rac 
for  the  random  fluctuations,  then  the  theory  given  in  I  may 
not  apply. 

We  conclude  in  this  section  that  it  is  reasonable  to  apply 
the  one-dimensional  fluid  theory  reviewed  in  Sec.  II,  with  the 
qualifications  stated  above,  to  turbulent  ion  heating  in  the 
long-range  potential  region  away  from  the  edges  of  the 
small-scale,  quasistatic  potential  structures  that  occur  in  the 
downward  Birkeland  current  sheets. 

IV.  CHARACTERISTICS  OF  FAST  SATELLITE 
OBSERVATIONS  FOR  SOME  DOWNWARD  CURRENT 
REGIONS 

A  number  of  different  plasma  and  wave  processes  occur 
in  downward  current  regions.  We  have  selected  three  FAST 
data  sets  (satellite  passes  in  a  downward  current  region)  with 
similar  properties.  These  properties  constitute  the  assumed 
physical  model  which  we  discuss  in  this  section.  For  a  more 
general  discussion  of  downward  current  regions  (see  Refs.  5 
and  10). 

For  downward  Birkeland  current  regions  in  the  Earth’s 
auroral  zone,  FAST  satellite  measurements  show  that  there 
are  often  four  plasma  constituents  present.  There  are  large 
densities  of  low-  to  intermediate-energy  (<few  keV)  elec¬ 
trons  and  ions  with  a  drift  velocity  up  the  flux  tube. 
The  electrons  are  field-aligned  (Tell^>Te±)  and  thermalized 
(ue<vel u), 10  and  the  ions  have  a  “conic”  shape  in  v±-V\\ 
phase  space,12  where  Ti±  >  Ti{l.  There  is  often  a  small  density 


Downloaded  21  Nov  2006  to  146.153.144.35.  Redistribution  subject  to  AIP  license  or  copyright,  see  http://pop.aip.org/pop/copyright.jsp 


1 1 2902-6  Jasperse  et  al. 


Phys.  Plasmas  13,  112902  (2006) 


of  high-energy  (5-30  keV)  ions  moving  down  the  flux  tube. 
Sometimes,  there  is  a  small  density  of  high-energy  (>few 
keV)  electrons  also  moving  down  the  flux  tube.  Since  the 
densities  of  the  high-energy  particles  at  FAST  satellite  alti¬ 
tudes  are  small,10  and  since  nearly  all  the  current  is  carried 
by  the  low-  to  intermediate-energy  electrons,  we  neglect  the 
high-energy  particles  and  consider  a  two-constituent  plasma 
consisting  of  the  low-  to  intermediate-energy  electrons  and 
ions. 

Satellite  data,  particle  simulations,  and  some  analytic  re¬ 
sults  show  that  there  are  at  least  four  main  altitude  regions 
for  downward  current  sheets:  (1)  a  region  below  the  double 
layer  (DL)  region  where  a  downgoing  ion  beam  has  been 
observed;30  (2)  one  or  more  strong  laminar  double  layers 
each  with  a  width  of  about  10\Df,-20\De,  at  the  bottom  of  the 
acceleration  region;13  31  (3)  a  transition  region  (TR)  just 
above  the  DL  region  with  a  width  of  about  \00\^e-200\De 
where  the  upgoing  electron  beam  is  thermalized;13,19  and  (4) 
a  long-range  potential  region  (LRPR)  that  extends  from  the 
top  of  the  TR  to  several  Earth  radii  and  beyond.12,14’23  24 
Statistically,  the  bottom  of  the  acceleration  region  (the  DL 
region)  occurs  during  the  winter  months  from  about  1600  to 
2500  km.32  The  LRPR  extends  from  the  top  of  the  TR  to  at 
least  4.3 RE  (geocentric)14  and  probably  beyond.  CLUSTER 
data14  also  show  that  the  downward  current  sheets  remain 
quasi  steady  on  time  scales  on  the  order  of  a  few  hundred 
seconds.  At  FAST  altitudes,  observations  of  the  VLF  saucers, 
sometimes  a  signature  of  the  downward  current  sheets,  indi¬ 
cate  that  the  downward  sheets  remain  quasisteady  on  time 
scales  greater  than  several  tens  of  seconds.10  The  turbulence 
in  the  TR  is  often  strong;19  in  the  LRPR  it  is  often  weak. 

FAST  observations  of  the  fluctuating  electric  field  in  the 
LRPR  show  that  there  are  at  least  two  kinds  of  turbulence 
associated  with  downward  auroral-current  regions:5  broad¬ 
band  extremely  low  frequency  (BBELF)12  and  electron  soli¬ 
tary  wave  (ESW).11  In  this  paper,  we  consider  only  the 
former.  In  a  brief  examination  of  FAST  data,  we  found  three 
satellite  passes  where  ion  heating  by  ESW  turbulence  is  neg¬ 
ligible,  but  strong  BBELF  fluctuations  are  present.  These 
three  cases  show  a  frequency  spectrum  of  electric-field  fluc¬ 
tuations  in  the  ELF  range  in  the  spacecraft  frame  of  refer¬ 
ence.  The  spectrum  consists  of  two  components:  a  weak, 
roughly  power-law  background  and  broad  structured  peaks 
in  the  spectral  density  ordered  by  the  proton  cyclotron  fre¬ 
quency.  The  peaks  are  due  to  temporal  fluctuations  of  the 
electric  field  in  the  ion  frame,  while  the  power-law  back¬ 
ground  is  due  to  Doppler-shifted,  spatially  irregular,  electric- 
field  structures  that  are  stationary  in  the  ion  frame.33  Sub¬ 
tracting  the  two  signals,  we  are  left  with  what  we  interpret  as 
real,  time-dependent  turbulence  in  the  ion  frame.34  In  Fig.  1, 
we  show  FAST  satellite  data  at  UT=  14:37:45.0274  for  the 
spectral  density  of  the  longitudinal  electric-field  fluctuations 
for  the  downward  current  sheet  beginning  at 
UT=1997-01-23/14:37:40  and  extending  for  about  10  s.  In 
panel  A,  we  show  the  satellite  observed  spectrum  and  in 
panel  B,  the  subtracted  spectrum.  The  vertical  dashed  lines 
denote  the  proton  gyrofrequency  and  its  first  two  harmonics 
in  Hertz.  In  each  of  the  three  data  sets  that  we  have  studied, 


f(Hz) 


FIG.  1.  FAST  observations  of  the  spectral  density  of  the  longitudinal 
electric-field  fluctuations  at  UT  1997-01-23/14:37:45.0274:  (A)  observed 
spectrum  and  (B)  subtracted  spectrum,  which  is  interpreted  as  time- 
dependent  electrostatic  turbulence.  Vertical  dashed  lines  denote  the  proton 
gyrofrequency  and  its  first  two  harmonics.  The  crossed-dashed  lines  in  panel 
B  show  the  “symmetric  window  function”  approximation  discussed  in 
Sec.  VI. 


we  have  determined  that  the  turbulence  is  electrostatic,  as 
discussed  in  Sec.  III. 

From  Eq.  (20)  of  I, 

w*  +Wel=-\  Wa  +  W,J.  +  jJ-[(SE2±)  +  (SE~)] 

dt  0  7T 

(24) 

FAST  satellite  data  in  the  LRPR  show  that  Wel  =  0  and  that 
Wi±  >  Wj^  0.  Since  we  are  considering  a  situation  where 
the  ions  gain  energy  from  the  turbulent  wave-ion  interaction, 
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and  since  the  turbulence  is  generated  by  a  plasma  instability, 
then  for  early  times 


.  d  1  dB 

rtiuX  — Ti\  - - r.i  =  WiL. 

'  '\ds  '±  Bds  ,LI  ,x 


(28) 


rr-[<®i>+<«ii>]>o.  (M) 

at  0  7T 

It  follows  from  (24)  and  (25)  that  Well  is  negative.  The  elec¬ 
trons  lose  energy  as  a  result  of  their  parallel  motion,  partly  to 
the  heating  of  the  ions  and  partly  to  the  unstable  growing 
electric-field  fluctuations.  If  a  quasisteady  solution  exists, 
then  the  wave  term  in  (24)  is  zero  and  IVd|  =  -(\V/1  +  VV,n). 
Since  WiL  +  Wi{l  is  positive,  then  in  the  quasisteady  state,  We[l 
is  also  negative. 

Therefore,  the  physical  model  that  we  adopt  for  the 
LRPR  is  that  of  a  magnetic-field-aligned,  U-shaped,  quasi¬ 
static  potential  structure  in  which  the  potential  slowly 
changes  as  the  distance  along  B  changes.5 10  14  We  consider  a 
flux  tube  in  the  U-shaped  region  away  from  any  small-scale, 
electrostatic  shocks  so  that  the  E±  X  B  drift  velocity  is  small 
compared  to  the  proton  drift  velocity  along  B,  as  discussed  in 
Sec.  III.  For  data  sets  of  the  above  type,  we  argue  that  it  is 
the  electrostatic  turbulence  that  heats  the  protons  by  cyclo¬ 
tron  resonance  perpendicular  to  B  in  the  LRPR. 

In  the  rest  of  this  paper,  we  assume  that  a  solution  of 
(l)-(5)  and  (8)— ( 1 0)  exists  that  is  quasisteady  on  a  time  scale 
of  a  few  hundred  seconds. 


V.  CLOSURE  APPROXIMATION 

FOR  THE  PERPENDICULAR  ION  DYNAMICS 

In  this  section,  we  discuss  the  issue  of  closure.  In  order 
to  solve  for  the  perpendicular  ion-temperature  profile  Ti±(s), 
we  must  specify  a  closure  condition  for  the  perpendicular  ion 
dynamics.  In  this  paper,  we  choose  the  closure  qi±(s) 
=  uJis)Ti±(s)  [see  (29)  below].  In  this  section,  we  also  dis¬ 
cuss  Chew,  Goldberger,  and  Low35  (CGL)  closure,  where  it 
is  assumed  that  Ti±(s)  ^B(s),  and  show  that  this  assumption 
is  not  valid  when  turbulence  occurs  on  downward  Birkeland 
current  sheets.  For  a  general  discussion  of  CGL  closure  (see 
Ref.  36). 

Consider  the  perpendicular  ion  dynamics.  Equation  (4) 
for  the  upward  drifting  ions  when  there  are  no  fluctuations 
present  is 

jniTil  +  B2y(niqil/B2)  =  0.  (26) 

dt  ds 

If  we  choose  the  closure  q{  then  for  low-frequency 

or  for  quasisteady  phenomena,  8TiLl  dt  is  negligible,  and  we 
would  have 


1  dB 
B  ds  7,1 


=  0, 


(27) 


where  we  have  used  (1).  This  closure  gives  the  CGL  result 
(Ti±~B)  when  the  fluctuations  are  negligible.  However, 
when  turbulence  is  present,  this  closure  gives  the  equation 


For  Uj>0  and  for  Wil/niui>-B~l(dB/ds)Til> 0,  the  ions 
will  be  heated  perpendicular  to  B  as  they  move  up  the 
B-field  line.  Therefore,  we  argue  that  it  is  reasonable  to 
choose  the  above  closure  relation  when  turbulent  ion  heating 
perpendicular  to  B  is  present.  Thus, 

<7;±  =  UiTi± .  (29) 


We  consider  an  electron  single-ion  plasma  with  closure  given 
by  (29)  and  approximate  (28)  by  choosing  an  average  value 
for  uh  denoted  by  uh  over  the  interval  jj 


\_dB 
B  ds 


= wa¬ 


rn 


Here,  wi±  is  the  perpendicular  ion-heating  rate  per  particle 
denoted  by  wil=nJlWi±. 

We  wish  to  use  the  results  of  Ref.  37  for  the  ion- 
temperature  relaxation  rate  to  estimate  the  perpendicular  ion¬ 
cooling  rate  for  a  nonturbulent  (collisional)  electron-proton 
plasma  with  an  ion  temperature  anisotropy  (Ti±>  T^)  in  the 
limit  as  B— ►O.  The  perpendicular  ion-temperature  relaxation 
rate  l/7?±  is  defined  as 

dT,L/dt  =  -(TiL-TJ/^x.  (31) 

From  Ref.  37,  we  find  that  at  FAST  altitudes 

1/t?x~3.1  X  10_,°  s_1 .  (32) 

Since  Ti±>TiihTi±  is  cooled  by  the  thermal  fluctuations  and 
the  perpendicular  ion-cooling  rate  per  particle  is  approxi¬ 
mately 


wi±{qu)  ~  -  3.1  X  1(T8  eV/s.  (33) 

Using  Ref.  37,  we  have  also  estimated  the  electron-ion  re¬ 
laxation  rate  for  quiescent  conditions  and  found  it  was  small 
compared  to  1/t^,  so  we  neglect  it. 

We  now  wish  to  compare  the  quiescent  cooling  rate 
per  particle  given  by  (33)  to  the  adiabatic  cooling  rate  per 
particle  given  by  the  second  term  on  the  left-hand  side 
of  (30).  For  FAST  satellite  altitudes,  B(s)  =  B0(a/s)2, 
0==6.4X  103  km,  5=  1.04 X  104  km,  i7;  =  50  km/s,  and 
fi±  =  150  eV,  so  we  have 

1  dB- 

w,“— T|X  =  -  2.2  eV/s.  (34) 

B  ds 

We  conclude  that  \wi±(qu)\  is  indeed  very  small  compared 

to  \uiB~](dB/ds)Til  \  and  that  fluctuations  at  the  thermal  level 
play  no  role  in  cooling  (or  heating)  the  ions  in  the  Earth’s 
magnetosphere.  The  above  argument  provides  a  justification 
that  the  CGL  closure  Ta  1  ~  B  is  valid  for  quiescent  (nontur¬ 
bulent)  magnetospheric  plasma  in  a  nonuniform  B  field. 
However,  when  there  is  turbulent  ion  heating  present  in  the 
plasma,  the  CGL  closure  Ti±  ~  B  is  not  valid. 
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VI.  ESTIMATE  OF  THE  TURBULENT  PERPENDICULAR 
ION-HEATING  RATE 

In  the  last  section,  we  estimated  the  adiabatic,  perpen¬ 
dicular  ion-cooling  rate  per  particle  for  downward  current 
sheets  in  the  LRPR  and  found  that  it  was  about  -2.2  eV/s  at 
FAST  altitudes.  In  this  section,  we  will  find  an  approximate 
formula  for  the  turbulent,  perpendicular  ion-heating  rate  per 
unit  volume  for  quasi-steady  conditions  using  (10).  It  is 
given  by  (49)  below  and  contains  what  we  call  the  “finite  ion 
gyroradius  effect,”  which  acts  to  reduce  the  perpendicular 
ion-heating  rate  for  large  values  of  s. 

From  ( 1 0)  using  symmetry,  we  see  that  Wi±  is  given  by 
a  four-dimensional  integral  and  an  infinite  sum  on  n.  If  we 
had  satellite  data  or  theoretical  formulas  for  (|££2|(.r,k,a>)) 
and  e(s;k,<o)  for. the  turbulent  plasma  for  quasisteady  con¬ 
ditions,  we  could  use  those  expressions  together  with  satel¬ 
lite  measurements  of  fi(syv±  ,yt|)  to  evaluate  W/x($).  How¬ 
ever,  since  such  detailed  information  is  lacking,  we  need  to 
make  a  number  of  assumptions  and  approximations  in  order 
to  estimate  VV/1(5).  These  assumptions  and  approximations 
are  discussed  in  this  section. 

When  an  instability  is  present  in  the  plasma,  it  is  pos¬ 
sible  for  the  turbulent  fluctuations  to  grow  to  be  many  orders 
of  magnitude  larger  than  those  of  the  thermal  fluctuations. 
Since  satellite  measurements  show  electrostatic  turbulence  in 
downward  current  regions,  we  may  infer  that  one  or  more 
instabilities  operate  on  the  flux  tube  and  are  sustained  by 
maintaining  contact  with  some  source  of  free  energy.  The 
expression  for  Wi±  is  defined  in  (6).  Using  (B24)  of  I,  which 


relates  Ffa  to  by  the  formula  F{=(1/2)(<?/<?v)-I>£,  we  see 
that  the  expression  for  Wi±  given  by  (6),  (7),  and  (10)  may 
be  rewritten  as 


Equation  (35)  has  the  form  of  a  friction  piece  given  by  the 
first  term  in  the  brackets  plus  a  diffusion  piece  given  by  the 
second  term.  It  can  be  shown  that  the  first  term  corresponds 
to  the  spontaneous  emission  of  plasma  waves  by  the  ions  and 
makes  a  negative  contribution  to  Wily  and  the  second  term 
corresponds  to  the  stimulated  absorption  and  emission  of 
plasma  waves  by  the  ions  and  makes  a  positive  contribution 
to  Wi±.  We  now  introduce  the  “diffusion  approximation  for 
the  perpendicular  ion  dynamics,”  where  we  assume  that  the 
perpendicular  ion  energy  loss  by  spontaneous  emission  is 
small  compared  to  the  perpendicular  ion  energy  gain  by 
stimulated  absorption;  therefore  we  neglect  the  first  term  in 
(35)  compared  to  the  second.  This  approximation  represents 
an  upper  bound  on  the  perpendicular  ion-heating  rate.  Since 
the  ion  velocities  are  small  compared  to  the  electron  veloci¬ 
ties,  we  argue  that  this  approximation  should  give  a  tight 
upper  bound  for  Wi±;  i.e.,  one  that  is  close  to  the  exact  value. 
From  (10),  we  therefore  neglect  the  term  proportional  to 
Im(l/e)  to  obtain  for  quasisteady  conditions 


Wi±(s) 


=  J  d3vft(s,v L,v^  J  J 


(36) 


FAST  satellite  data  for  the  ion-distribution  function  show  that  it  has  a  conic  shape  in  space.  Therefore,  we  are  led 

to  approximate  /,  by  a  “pancake-shaped  function” 


fi(s,V±,Vu)  =fi 


syv±y 


fall  -  «/)2 

2»i  J’ 


(37) 


where  /,  is  a  spread-out  function  of  vL{TiL  >  7^-n)  and  a  peaked  function  around  v  {=uh  where  w,  is  positive  for  ions  flowing  up 
the  flux  tube.  This  is  a  reasonable  approximation  for  /,  at  lower  altitudes  before  it  becomes  folded  up  by  the  action  of  the 
mirror  force  term  in  (2)  at  higher  altitudes.  We  now  integrate  (36)  on  v  \  and  use  integration  by  parts  on  o)  in  the  second  term 
of  (36)  to  obtain 


(<o  -  n(lj  -  knUj)2 


(g>  -  nil,  -  kjUj)2 
2kjvl 


2k- 


x«fa;_i(^,)  -•/;+,(£)]  -/'  s,v±, 


2klvl 


)(w-nlli-*l|Ml)^fa()  [<|<S£2|(s,k,w)).  (38) 

k\\Vi\\ 
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where  /'  is  the  derivative  of/,  with  respect  to  its  third  argu¬ 
ment.  The  following  ordering  holds  for  downward  currents: 


0  <  ui  <  Vi\\  and  0  <  |&J  <  k±y 


(39) 


fi(s,v1,vl)  =  ni(s) 


1 

2lw?± 


S, 


(i'll  ~  “.)2 

2vf„ 


(40) 


where  and  k±  are  average  or  “characteristic”  wavenumbers 
for  the  turbulence.  Using  the  approximation  that  /,  is  a 
peaked  function  of  v\\  around  w„  the  ordering  in  (39),  and 
k±vil/Cli<  1,  we  see  that  for  Ti±/T^  not  too  large,  the  sec¬ 
ond  term  in  (38)  is  small  compared  to  the  first.  Therefore,  we 
may  neglect  it. 

In  order  to  make  further  progress  with  the  estimation  of 
Wi±,  we  need  to  make  another  approximation  for  the  form  of 
/,.  For  our  purposes  in  this  section,  we  may  approximate/,  in 
(36)  by  the  product  of  a  Maxwellian  in  v±  times  the  peaked 
function  of  yy: 


sr 


fan  -  “,)2 


2v 


ill  J 


J 


drk  y 

Mix 


XAB(ft){<|«2|(j;k,«ni  +  *||U||)>  +  (|<5£2|(5  ;k,- 


where  viL=vi±(s),  and  7,  {.v,(i>  —u,)2/2vfy  is  normalized  to 
1 .  As  we  will  see  below,  the  main  effect  of  this  approxima¬ 
tion  is  to  cut  off  the  v± -integration  and  introduce  the  func¬ 
tion  (2 /#■) A !(/?,-)  into  the  formula  for  vV,±(j).  Here, 
and  A„(jc)=exp(-jc)/n(x),  where  ln(x)  is 
the  Bessel  function  of  imaginary  argument  of  order  n.  A 
more  accurate  treatment  of  /,,  for  example,  a  two- 
dimensional  orthonormal  expansion  in  the  variables  v±  and 
(yn-M/)2/2yJ|,  would  give  a  somewhat  different  cutoff  for  the 
v± -integration,  but  should  not  change  the  main  result  in  this 
section.  Substituting  (40)  into  (36)  and  evaluating  the  v±- 
and  co-integrations,  we  see  that  the  first  term  of  (36)  is 


fifl, +  *„!;„)>}.  (41) 


r 


From  Sec.  II,  we  see  that  (|££2|(5;-k,-co)) 
=  (|££2|Cs;k,co)).  This  means  that  there  is  as  much  spectral 
density  at  -k  and  -co  as  there  is  at  k  and  co.  As  we  mentioned 
earlier,  if  we  had  a  satellite  measurement  of  (\SE2\  ($;k,co)), 
we  could  then  use  it  to  evaluate  the  ^-integration  in  (41). 
Since  we  are,  at  present,  lacking  such  detailed  information 
about  the  spectral  density,  and  in  order  to  make  further 
progress  with  the  approximate  evaluation  of  (41),  we  assume 
that  the  k  and  co  variables  in  the  spectral  density  are  approxi¬ 
mately  separable.  In  general,  we  do  not  expect  the  spectral 
density  to  be  exactly  separable  on  k  and  to,  even  though  the 
plasma  may  be  very  turbulent  and  have  no  linear  dispersion 
relation.  We  argue  that  this  approximation  should  be  good 
enough  for  the  approximate  analysis  in  this  section.  There¬ 
fore,  we  make  the  “separation  of  k  and  co  variables”  approxi¬ 
mation  and  write 

(\SE2\(s;k,<o))  =  <|^2|(js;k))<|^E2|(^;o>)>,  (42) 

where  (|<5£2|(s;k))  is  normalized  to  1.  This  approximation 
means  that  (|<5£2|(s;-k)}s(|<5£2|(,s;k))  and  (|<5£2|(j;-w)) 
=  (|<5£2|(s;w)>. 

From  satellite  data  for  downward  currents,  we  also 
find  that  the  values  of  (|<5E2|($;cu))  as  a  function  of  (o  are 
much  smaller  near  successive  values  of  (i)=±n£lh  where 


n- 2,3,...,  than  the  values  of  (|<$£2|($;cu))  near  cu=±fl,.34 
For  example,  see  Fig.  1.  We  are  led  to  assume  that 
(|<$E2|Cs;o)))  decreases  sufficiently  fast  for  large  |a>|  so  that 
the  low-order,  co-moments  of  (|<5E2|(s;co))  exist.  As  a  result, 
we  are  also  led  to  make  the  “symmetric  window  function” 
approximation  for  (|££2|(s;co)),  which  is  defined  as 

[0,  otherwise, 

(43) 

where 

/•+  00 

(|<5£2|(s))a„  =  (2a>0)"'  I  da>(|<S£2|(.r;a>)).  (44) 

J  — X 

The  meaning  of  this  approximation  is  that  the  value  of 
(|<$E2|(j;co))  as  a  function  of  co  is  replaced  by  its  average 
(constant)  value  for  -co0<co<  +co0  and  by  zero  for 
|co|>co0,  where  co0  is  the  half-width  of  (|<$E2|(.y;co))  in  co 
space.  The  half-width  in  co  space  is  defined  as 


Downloaded  21  Nov  2006  to  146.153.144.35.  Redistribution  subject  to  AIP  license  or  copyright,  see  http://pop.aip.org/pop/copyright.jsp 


1 1 2902- 1 0  Jasperse  et  al. 


Phys.  Plasmas  13,  112902  (2006) 


co0=  I  d(i)(o(\SE2\(s\(x)))/  I  d(j){\SE2\(s\(o)).  (45) 

Jo  Jo 

The  symmetric  window  function  as  a  function  of  the  fre¬ 
quency  /(=cu/27t)  is  shown  graphically  in  Fig.  1  by  the 
crossed-dashed  lines  in  panel  B,  where  fo=o>0/ 2tt=257  Hz 
and  (\SE2\)av=2A  X  10"7  (V/m)2  Hz-1  at  the  satellite  alti¬ 
tude.  Here,  for  the  three  FAST  satellite  data  sets  that  we  have 
studied,  we  find  that  co0  is  a  number  that  varies  from  about 
1.1ft,-  to  1.4ft/.  From  satellite  data  for  the  ion  velocity  mo¬ 
ments  and  from  estimates  of  typical  values  for  |fcj,  the  fol¬ 
lowing  ordering  of  parameters  for  downward  current  regions 
at  FAST  altitudes  holds:  0 <ui<vi]l  and  l^n | l?/h < ft,  — ' co0.  Us¬ 
ing  this  ordering  of  parameters  and  the  symmetric  window 
function  approximation,  we  see  that  for  the  data  sets  consid¬ 
ered,  the  sum  in  (41)  is  dominated  by  the  n=  1  term  and, 
since  7/h{5,(u||-M/)2/2i;^|}  is  normalized  to  one,  the  integra¬ 
tion  on  v\\  gives 


X  (  ^- )  A !  (j8,)(|  <SE2|  (* ;  k)> .  (46) 

We  now  assume  that  (|<5E2|(j;k))  decreases  sufficiently 
fast  for  large  kL  and  |&n|,  so  that  its  low-order  klt  and  kl{ 
moments  exist.  Therefore,  we  are  also  led  to  approximate  the 
k±-  and  ^-integrations  by  making  “half-width-in-fc-space” 
approximations.  This  gives 


since  (|<5E2|(.s;k))  has  been  normalized  to  1.  Here, 
Pi0=k2± qV^/CI2  and  fco=*±o+^o-  The  half-widths  in  kl{  and 
k±  space  are  defined  as 


37.  In  this  case,  the  integrations  were  cut  off  using  the 
“dominant  term”  approximation  of  Chandrasekhar.^  If  the 
low-order  k±  and  k\\  moments  and  the  normalization  of 
(|<$E2|(s;k))  do  not  exist,  then  some  other  appropriate  way 
of  estimating  “characteristic”  values  for  k±  and  k{l  needs  to 
be  found,  so  that  k±0  and  /c|10  can  be  estimated.  We  also  note 
that  in  (49),  k±0  and  km  are  functions  of  s,  and  that  vi±  and 
ft/,  which  appear  in  /3l0,  are  also  functions  of  s. 

In  making  the  above  assumptions  and  approximations, 
we  have  reduced  the  four-dimensional  integrations  and  the 
infinite  sum  in  (10)  to  a  simple  formula.  In  (49),  we  see  that 
knowledge  of  the  average  spectral  density  in  w-space, 
(|££2|(s))fly,  the  half- width  in  a>-space,  co0(s),  the  half-widths 
in  k±-  and  kr space,  k10(s),  and  km(s),  ft,(j),  and  vi±(s)  de¬ 
fine  VT/j^s).  In  (49),  the  function  (2//3<0)A is  a  mono- 
tonically  decreasing  function  of  its  argument/*  /3i0 
=k2LOv2±/Cl2.  We  note  that,  as  the  ions  are  heated  perpen¬ 
dicular  to  B  and  move  up  the  flux  tube,  ui±  increases  from  a 
small  value  at  the  bottom  of  the  heating  region  where 
k±0vi±/Elj<\,  to  a  larger  value  at  high  altitudes  where 
^xo^ij./ ftf >  1.  In  this  way,  the  “finite  ion  gyroradius  effect” 
will  act  to  reduce  the  perpendicular  ion-heating  rate  for  suf¬ 
ficiently  large  s.  From  FAST  satellite  data  near  4000  km 
from  the  Earth’s  surface,  we  estimate  that  k2±o>k20  and  (3i0 
<1.  In  this  limit,3*  we  find  that  (2/j8/0)A1(j8/0)sl.  There¬ 
fore,  a  rough  estimate  for  (49)  at  FAST  and  lower  altitudes  is 

(50) 

Here,  the  units  of  (|<S£2|(.s))fli;  have  been  changed  from 
(V/m)2  (angular  frequency)-1  to  (V/m)2  Hz-1.  We  note  that 
(50)  is  1/2  the  value  obtained  heuristically21  from  dimen¬ 
sional  arguments  for  the  perpendicular  ion-heating  rate  per 
particle  for  electromagnetic  ion-cyclotron  turbulence  in 
downward  auroral-current  regions.  FAST  satellite  measure¬ 
ments  show  that  for  the  downward  current  sheet  beginning  at 
UT  1997-01-23/14:37:40  and  extending  for  about  10  s, 
(\SE2\)av~~2A  X  10-7  (V/m)2  Hz-1,  so  we  obtain 


*ii 


(|<$E2|(.s;&i,*il)>. 


(48) 


Using  these  results,  we  obtain 

WiAs)  s  )(^)A,(A°)<|^2|w>a„. 

(49) 

Since  {k\  / 2A:2) ( 2 / /5#) A  j  (>S/)  in  (46)  is  a  bounded  func¬ 
tion  of  k±  and  k\\,  and  (\SE2\(s\k1,kll))  is  a  positive-valued 
function  of  k±  and  klly  this  approximation  makes  sense  in 
light  of  the  first  mean  value  theorem  for  infinite,  one¬ 
dimensional  integrals  (see  Ref.  38,  Sec.  12),  provided  that 
there  is  no  singularity  in  the  iterated  kL-  and  ^-integrations. 
We  note  that  the  /.'-integrations  in  (46)-(48)  could  contain  a 
logarithmic  singularity  for  large  k,  as  was  found  in  the  ex¬ 
pressions  for  the  relaxation  rates  for  quiescent  plasma  in  Ref. 


wi±ao  ~  5.0  eV/s.  (51) 

We  see  from  (34)  that  wiLav>\uiB“x(dB/ds)fiJL\  so  that 
Ti±(s)  will  increase  as  the  ions  move  up  the  flux  tube.  How¬ 
ever,  as  shown  by  (49),  this  increase  will  not  continue  indefi¬ 
nitely,  and  as  the  ions  heat  up  and  k±QVil/Cli^  1,  Wi±(s) 
will  begin  to  decrease  with  increasing  s  by  the  finite  ion 
gyroradius  effect. 

A  comparison  of  (33)  for  the  nonturbulent,  very  weakly 
collisional  plasma  and  (51)  for  turbulent  plasma  shows  that 
the  electric-field  fluctuations  have  grown  by  about  eight 
orders  of  magnitude  due  to  the  presence  of  one  or  more 
instabilities  and  a  source  of  free  energy  on  the  downward, 
Birkeland  current  sheet. 

VII.  CALCULATION  OF  THE  PERPENDICULAR 
ION-TEMPERATURE  PROFILE  AT  LOW  ALTITUDES 

In  this  section,  we  solve  for  the  low-altitude,  perpendicu¬ 
lar  ion  temperature  profile,  Ti±(s)  in  terms  of  wilav(s)  and 
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w,($)  for  quasisteady  conditions.  The  formula  is  given  by  (56) 
below  and  parametrized  calculations  of  7(1($)  are  also  given 
in  Fig.  2. 

As  discussed  in  Sec.  V,  the  quasisteady  equation  for 
Ti±(s)  from  (28)  is 

dTiL(s)lds  -  {B(s)~'dB(s)/ds}Ti±(s)  =  wi±av(s)/Ui(s), 

(52) 


Wu«,(J)  =  nj(s)  1  VV/1ot(j)  . 


(53) 


Equation  (52)  may  be  solved  by  the  method  of  integrating 
factors.  The  general  solution  is 


Ti±(s)  =  ex  p 


B  uB/ds)  , 


xjds  wi±av(s)/Ui{s) 

Xexpj+J  ds(-  +  C 

Xexpj-  J  ds(-  j. 


(54) 


where  C  is  a  constant  that  specifies  a  boundary  value  for 
Ti±(s).  We  may  verify  that  (54)  is  a  solution  of  (52)  by 
substitution.  From  our  study  in  Ref.  23,  Fig.  1,  we  found  that 
(|<5£2|($))fli;  increases  with  5  (as  ft,-  decreases)  and  is  ap¬ 
proximately  proportional  to  the  increase  in  m,($)  with  $,  so 
that  wilav(s)/ ufa)  was  nearly  constant  for  relatively  small 
distances,  $2-$i  < few  X  103  km.  At  these  altitudes, 
kl0vil /ft,-<  1,  so  that  the  heating  rate  was  not  reduced  by 
the  finite  ion  gyroradius  effect.  Therefore,  we  make  the  ap¬ 
proximation 


=  [wafl®(JiVM/(Ji)](^i)y»  (55) 

where  -0.5^  y ^  +0.5.  In  addition,  - B~]dB/ds  =  3/s ,  so  the 
solution  (54)  is 


Ti±(s)  =  (si/s)3TiL(s\)  +  {Wi±av(s\)s~\  y/Ui(si)(4  +  y)} 

X{sUy-s\+y(Sl/s)3},  (56) 


where  Wilav(s\)  =  {q]l^m^(\8E2\(sx))av. 

In  order  to  show  that  (56)  gives  the  correct  result  for 
Ti±(s),  we  would  need  measurements  by  two  satellites  at 
different  altitudes  on  the  same  flux  tube  at  the  same  time. 
One  satellite  would  give  the  boundary  condition  for  Ti±  at 
(the  bottom  of  the  LRPR)  and  values  for  (\SE2\)av  and  w„ 
and  the  other  satellite  would  give  the  heated  values  of  Ti±  at 
the  satellite  altitude  s2.  Since  such  coincident  measurements 
are  not  available,  in  this  section  we  present  parametrized 
calculations  of  Tu  as  a  function  of  vvllflu  and  m,  at  s2  and  the 
boundary  condition  for  Ti±  at  $j.  Satellite  data  for  typical 
downward  current  sheets  show  that  the  ions  are  heated  from 
about  20  to  80  eV  at  the  bottom  of  the  LRPR  to  about  100  to 
200  eV  at  the  satellite  altitude  near  4000  km.5  Satellite  data 
also  show  that  the  bottom  of  the  acceleration  region  for  win¬ 
ter  months  varies  from  about  1600  to  2500  km/2  Therefore, 


FIG.  2.  Solutions  of  (56)  for  Ti±(sf)  in  eV  as  a  function  of  s'  from  2000  to 
6400  km  for  7'/1(jJ)=50  eV,  y=0,  Uj(s£=50  km/s,  for  four  values  of 
vv,lfly(jJ)  given  by  50,  5.0,  0.5,  and  0  in  eV/s.  The  vertical  dashed  line 
denotes  the  satellite  altitude  at  $2=4130  km. 


we  choose  ^  =  8400  km  ($J  =  2000),  Ti±(s[)=50  eV,  and  s2 
=  10  530  km  ($2=4130)  as  the  satellite  altitude.  Here,  $'  is 
the  distance  along  B  from  the  surface  of  the  Earth,  $'=$-a, 
where  a  is  the  Earth’s  radius.  From  FAST  data,  we  make 
estimates  of  (\SE2\)av  and  w,  at  the  satellite  altitude,  s2,  and 
calculate  the  corresponding  value  of  wilav(s2^u^s2^  us^ng 
(50).  The  values  of  vv,ifli;($|)/w,($|)  are  inferred  from 
vv,lfl„($2)/w/(-S2)  by  using  (55)  for  a  choice  of  the  value  of  y. 
They  are  then  used  in  (56)  to  calculate  Ti±(s'). 

Calculations  of  Ti±(s')  in  eV  from  $'  =  2000  km  to  $' 
=  6400  km  are  shown  for  7’,±($J)  =  50  eV,  Uj(s2)=50  km/s, 
y=0,  and  for  four  values  of  Wi±av(s'2)  in  eV/s,  denoted  by 
the  dashed  and  solid  curves  in  Fig.  2.  The  vertical  dashed 
line  shows  the  satellite  altitude  at  $2=4130  km.  We  note  that 
when  wilaB(s^=0,  Tu(s')=Tu(s'1)[(s'1+a)/(s' +a)f 
~ /?($').  This  means  that  when  there  is  no  turbulence  present, 
the  CGL  closure  condition,  as  discussed  in  Section  V,  is  valid 
and  that  Tu  decreases  with  $',  as  shown  by  the  dashed  curve 
in  Fig.  2.  However,  when  turbulence  is  present,  Fig.  2  shows 
that  T,±  increases  with  $'  and  the  CGL  closure  condition  is 
no  longer  valid.  Note  that  the  curve  for  wilav  =  5  eV/s  pro¬ 
duces  a  Ti±  of  182  eV  at  the  satellite  altitude  and  is  typical 
of  the  three  FAST  data  sets  that  we  have  studied.  Note  also 
that  if  larger  heating  rates  are  present,  Ti±(s2)  in  the  few  keV 
range  can  occur.  Calculations  were  made  for  7= -0.5  and 
7= +0.5  and  the  values  for  Ta  ($2)  were  shifted  up  and 
down  by  about  4%  from  those  given  in  Fig.  2.  Calculations 
were  also  made  for  T,i($J)=20  eV  and  80  eV,  and  for 
$J  =  1600km  and  2500  km,  and  curves  similar  to  those 
shown  in  Fig.  2  were  obtained. 

The  results  presented  in  Fig.  2  can  be  understood  in 
simple  terms.  Equation  (52)  is  basically  a  rate  equation  for 
Ti±(s').  Consider  the  curve  in  Fig.  2  for  Ti±(sf)  for 
Wi±av(s 2)  =  5  eV/s.  From  (34),  we  see  that  the  adiabatic  cool¬ 
ing  rate  is  about  -2.2  eV/s  and,  since  the  turbulent  heating 
rate  is  5  eV/s,  the  net  heating  rate  is  about  2.8  eV/s.  From 
FAST  data,  the  average  upward  drift  of  the  ions  was 
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50  km/s,  so  it  takes  about  43  s  for  an  ion  to  travel  from 
sj  =  2000  km  to  $2=4130  km.  Thus,  the  average  energy  gain 
for  an  ion  is  43  sX2.8  eV/s=120  eV,  and  50  eV+120  eV 
=  170  eV.  This  compares  well  with  the  exact  value  of  182  eV 
obtained  from  (56)  for  y=0. 

VIII.  SUMMARY  AND  DISCUSSION 

The  main  results  of  this  paper  are  summarized  as  fol¬ 
lows. 

(1)  In  Sec.  II,  we  have  reviewed  a  new  fluid  theory  for  the 
Birkeland  current  system  in  the  guiding-center  and  gy- 
rotropic  approximation  derivable  from  the  ensemble- 
averaged  Vlasov-Maxwell  equations  that  includes  the 
effect  of  wave-particle  interactions  in  which  the  particles 
are  transported  in  one  spatial  dimension  (the  distance  s 
along  B),  but  the  turbulence  is  two-dimensional  since 
electric-field  fluctuations  that  are  both  perpendicular  and 
parallel  to  B  are  considered. 

(2)  When  the  turbulence  is  random  and  electrostatic,  we  ar¬ 
gue  that  the  one-dimensional  fluid  theory  may  be  ap¬ 
plied  in  the  long-range  potential  region  (LRPR)  of  a 
downward  Birkeland  current  sheet  on  flux  tubes  away 
from  the  edges  of  any  small-scale,  U-shaped,  quasistatic 
potential  structures,  so  that  the  Ei  X  B  drift  remains 
small  compared  to  the  ion  drift  along  B  (see  Sec.  III). 

(3)  In  Sec.  V,  we  use  qi±=uiTil  to  close  the  perpendicular 
ion  dynamics  and  present  an  argument  as  to  why  this  is 
a  reasonable  assumption. 

(4)  For  quasisteady  conditions,  we  have  reduced  the  expres¬ 
sion  for  Wi±  from  a  four-dimensional  integral  and  an 
infinite  sum  on  n ,  given  by  (10),  to  a  simple  formula, 
given  by  (49),  by  making  a  number  of  assumptions  and 
approximations.  They  are:  (a)  the  “diffusion  assumption 
for  the  perpendicular  ion  dynamics”;  (b)  the  approxima¬ 
tion  of  the  ion  conic  in  v±-V\\  space  by  a  “pancake¬ 
shaped  function”  with  the  vL  -dependence  given  by  a 
Maxwellian;  (c)  the  assumption  that  (|<$E2|($;k,a>))  is 
“separable  in  its  k  and  a)  variables”;  and  (d)  the  idea  that 
(|  <5E2 1  (,s;<i)))  can  be  approximated  by  a  “symmetric  win¬ 
dow  function”  and  that  average  or  “characteristic”  val¬ 
ues  for  kL  and  kl{  can  be  found  from  (|<5£2|($,k)).  We 
argue  that  (49)  gives  a  reasonable  upper  bound  for  WiL. 
Our  analysis  leads  to  the  idea  that  the  turbulent,  perpen¬ 
dicular  ion-heating  rate  does  not  increase  indefinitely 
with  increasing  s ,  but  is  reduced  for  large  5  by  the  “finite 
ion  gyroradius  effect”  (see  Sec.  VI). 

(5)  In  Sec.  VII,  we  gave  a  new  formula  for  the  low-altitude, 
perpendicular  ion-temperature  profile,  Ti±(s'),  and  also 
gave  parametrized  calculations  of  TiAs')  as  a  function 
of  Wi±m(s'l)  and  ui(s2 )  in  Fig-  2* 

There  are  other  mechanisms  by  which  magnetospheric 
ions  can  be  heated  perpendicular  to  B  on  auroral  field  lines 
beside  the  one  presented  in  this  paper.  For  a  discussion  of 
these  processes  (see  Ref.  5).  When  ESW  turbulence  is 
present  in  the  LRPR  of  downward  current  sheets,  it  is  well 
known  that  intense  ion  heating  perpendicular  to  B  occurs.11 
If  ion  heating  by  ESW  turbulence  is  a  random  process  and  if 


a  theory  for  the  spectral  density  of  the  longitudinal  electric- 
field  fluctuations  could  be  developed,  then  the  spectral  den¬ 
sity  could  be  used  to  find  wi±  for  ESW  turbulence,  and  cal¬ 
culations  of  TiL  as  a  function  of  s'  could  also  be  made. 

The  following  picture  for  the  downward  Birkeland  cur¬ 
rent  region  is  suggested  by  satellite  observations,5'14  particle 
simulations,16'20  and  the  analysis  in  this  paper  and  in  Refs.  1 
and  21-24.  As  the  current  density  increases  with  decreasing 


altitude,  one  or  more  strong  double  layers 


13,16-20,31 


develop 


at  an  altitude  that  ranges  during  winter  months  from  about 
1600  to  2500  km  from  the  Earth’s  surface.32  The  electrons 
are  accelerated  upward  by  the  downward  pointing,  self- 
consistent  Ey  field  and  some  of  the  ions  are  dragged  upward 
against  the  downward  pointing  EM  field  in  order  to  maintain 
quasineutrality.8’16'20  23  24  In  the  LRPR  when  ESW  turbu¬ 
lence  is  absent,  the  ions  are  heated  perpendicular  to  B  by 
cyclotron  resonance  with  electrostatic  turbulence  near  fly  and 
its  harmonics.1  The  ions  are  gaining  perpendicular  and  par¬ 
allel  energy  from  the  turbulence  and  the  electrons  are  losing 
parallel  energy  to  the  turbulence,  as  indicated  by  (24)  and 
(25)  of  Sec.  IV.  In  a  quasisteady  state,  the  amount  of  parallel 
energy  lost  by  the  electrons  is  equal  to  the  amount  of  per¬ 
pendicular  and  parallel  energy  gained  by  the  ions.1  In  the 
LRPR,  some  of  the  electrons  that  are  accelerated  upward  by 
Eh  are  driven  into  the  downward  direction  by  the  magnetized, 
turbulent,  electron-wave  interactions.1  Since  Wel  =  0,  the 
electrons  are  actually  cooled  adiabatically  in  the  perpendicu¬ 
lar  direction9’10  as  indicated  by  the  perpendicular  energy  bal¬ 
ance  equation  for  the  electrons.1’23’24  Since  W<>1=0,  these 
equations  for  the  electrons  are  similar  to  (26),  (27),  and  (29) 
of  this  paper.  In  the  LRPR,  this  leads  to  a  field-aligned 
(T«*TeJ  and  thermalized  (ue<ve\)  electron  distribution 
function  parallel  to  B, 1,10,23,24  and  a  heated  ion-distribution 
function  (7, ,  >  Tq)  in  the  directions  perpendicular  and  par¬ 
allel  to  B.1,10’12’23,  4  In  the  LRPR,  the  mirror  force  acts  on  the 
ions  [see  the  third  term  on  the  left-hand  side  of  (2)],  so  as  to 
fold  up  the  ion-distribution  function  into  a  conic-shaped  dis¬ 
tribution  in  v±-v\i  phase  space,  as  the  ions  move  up  the  flux 
tube.12  23,24  In  the  LRPR  for  quasisteady  conditions,  prelimi¬ 
nary  calculations  using  ( 1 )— (5)  and  (8)— ( 1 0)  show  that  the 
self-consistent  Ey  is  sustained  primarily  by  the  mirror  force 
(Alfven-Falthammar)  and  parallel  energy  divergence  terms 
in  the  momentum  balance  equation  given  by  (2), 1,1 2,23,24  and 


that  the  anomalous  resistivity,  which  is  proportional  to  Afdi, 
plays  a  minor  role.  The  above  scenario  is  consistent  with 
FAST  and  other  satellite  observations  for  downward  Birke¬ 
land  current  sheets.6"14’31’32 

We  recognize  that  a  two-dimensional  fluid  theory  in  the 
guiding-center  coordinate  system  needs  to  be  developed  so 
that  ExXB  drifts  near  the  edges  of  the  small  scale 
U-shaped,  quasistatic  potential  structures  may  be  taken  into 
account.  An  analysis  along  these  lines  is  currently  in 
progress  where  we  assume  that  the  geomagnetic  field  is  that 
of  a  dipole  and  use  dipolar  coordinates  to  express  the 
guiding-center  kinetic  equations,  Poisson’s  equation,  the  mo¬ 
ment  equations,  and  the  wave-particle  transfer  rates  for  the 
turbulent  plasma.  We  also  recognize  the  need  to  include  the 
effect  of  electromagnetic  turbulence. 
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It  is  the  main  conclusion  of  this  paper  that,  in  the  ab¬ 
sence  of  ESW  turbulence  when  structured  electrostatic  tur¬ 
bulence  near  and  its  harmonics  is  present  on  downward 
Birkeland  current  sheets  in  the  long-range  potential  region, 
the  ions  are  heated  perpendicular  to  B  by  cyclotron  reso¬ 
nance,  and  the  gyrotropic  guiding-center  fluid  theory  pre¬ 
sented  in  this  paper  gives  the  perpendicular  ion  temperature 
profile  at  low  altitudes  for  quasisteady  conditions. 
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